Dynamical behavior of interacting dark energy in loop quantum cosmology 
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The dynamical behaviors of interacting dark energy in loop quantum cosmology are discussed in 
this paper. Based on defining three dimensionless variables, we simplify the equations of the fixed 
points. The fixed points for interacting dark energy can be determined by the Friedmann equation 
coupled with the dynamical equations in Einstein cosmology. But in loop quantum cosmology, 
besides the Friedmann equation, the conversation equation also give a constrain on the fixed points. 
The difference of stability properties for the fixed points in loop quantum cosmology and the ones 
in Einstein cosmology also have been discussed. 
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I. INTRODUCTION 

Dark energy [lj is introduced to explain the expand- 
ing of our universe. We still have very little knowledge 
about it. One of the more interested candidate is cos- 
mological constant with the equation of state parame- 
ter u be = — 1< But the physical interpretation of the 
vacuum energy density is still faultiness. And the prob- 
lem that why the dark energy and dark matter are at 
the same order today is still unknown, this is so-called 
coincidence problem. To overcome this problem, many 
dynamical dark energy models have been introduced to 
replace the cosmological constant. But there are still 
many other approaches to explain the coincidence prob- 
lem. An interesting proposal is that the dark energy is 
interacting with dark matter. Many interacting models 
have been studied. The interacting term always is the 
function of the Hubble parameter H, the density of dark 
energy pde and the density of the dark matter pdm (for 
example, please see the second section of [2j and the ref- 
erences of it), or the function of T x , pde and pdm 0-Q- 

A general interaction between dark energy and dark 
matter is described by 

Pde + 3Hp DE {l +ojde) = Q, (1) 
Pdm + "AH p DM = -Q. (2) 

in which ujde is the equation of state parameter of dark 
energy. A constant lode and lode = —1 crossing were 
studied by many authors 0-111 • I n this paper we just 
consider lode is & negative constant. Q is the interacting 
term. And we will discuss two models: 



Model I Qi=3HCip D E 
Model II 



3HC2PDM, 



o2— n/~<Tj3— 2n ~n 

2 — o PDM: 



(3) 
(4) 



in which CijCvjC are some constants. Model I has been 
discussed by 0,0] in Einstein cosmology (EC), they con- 
sidered that C\ 7^ C2. But more special cases also have 



been studied, e.g., d = C 2 [lfj, d = [U|, C 2 = [T2J, 
in EC. For more simplicity, we will discuss C\ = C2 = C 
for the first model. The second model also has been dis- 
cussed in EC [H, [13], n is constant, and we will discuss 
n = 2. We will introduce three dimensionless variables to 
study those models in loop quantum cosmology (LQC), 
and discuss the difference and the common properties for 
dynamical system between LQC and in EC. 
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LQC [15Ml7l | is a canonical quantization of homoge- 
neous spactimes based upon techniques used in loop 
quantum gravity. Due to the homogeneous and isotropic 
spacetime, the phase space of LQC is simplifier then loop 
quantum gravity, e.g., the connection is determined by a 
single parameter called c and the triad is determined by 
p. The variables c and p are canonically conjugate with 
Poisson bracket {c,p} = ^p, in which 7 is the Barbero- 
Immirzi parameter. In LQC scenario, the initial singu- 
larity is instead of a bounce. So thanks for the quantum 
effect, the universe is initially contracting phase with the 
minimal but not zero volume, and then the quantum ef- 
fect drives it to the expanding phase. In the effecitve sce- 
nario, the behavior of LQC also can be well described, 
e.g., the big crunch instead of the big bounce [H, Il9j . 
People always consider two kinds of correction of effec- 
tive LQC, one is inverse volume correction, the other 
one is holonomy correction. In this paper, we just dis- 
cuss the holonomy correction. In this effective LQC, the 
Friedmann equation is modified (see Eqs. (f!M)) h It is easy 

to get the difference between the Friedmann equation in 

2 

EC and the one in LQC: it adds a — term to stan- 

Pc 

dard Friedmann equation which essentially encodes the 
discrete quantum geometric nature of spacetime [10, [HI ■ 
This means that the energy density will have a critical 
value, when the energy density is very nearly the critical 
density p c , the cosmology will have a bounce and then 
oscillates forever. Due to the modification of the Fried- 
mann equation, the dynamical behaviors of dark energy 
and dark matter are very different from the ones in EC. 
This factor was showed by many authors, for example, 
the phantom field (22|, the quintom field (23|, the in- 
teracting phantom field ;24| , the interacting dark energy 
with the interacting term Q = 3cH Pde Pdm , and 
the dark energy interacting with cold dark matter with- 



2 



out coupling to the baryonic matter 25]. In this paper, 
we will discuss the dynamical behavior of interacting dark 
energy. 

This paper is organized as following. In Sec. [TTl we 
will review the dynamical behavior for interacting dark 
energy in EC. And Sec. IIIII will show the dynamical be- 
havior of interacting dark energy in LQC. The numerical 
analysis will be showed in the Sec. IIVI And we get some 
conclusions in the last Sec. El For simplicity, we set 
8ttG=1. 



argued, if z = z{x, y), the above equations are closed and 
autonomous. The model I ^ and II ((H) both are in this 
case. 

For a autonomous system X' = /(X), the fixed points 
X* satisfy X' = 0. So, setting x' = 0,y' = of 
Eqs. (|llll2[) and simplifying them, it is easy to get the 
fixed points satisfy: 

x* + j/* = 1, (13a) 
— 3ujdex*(1 — a;*) + = 0, (13b) 



II. DYNAMICAL BEHAVIOR IN EC 



In EC, the Friedmann equation can be written as 
H 2 = U 



(5) 



in which p = pde + Pdm is the total energy density. 
Differentiating Eq.([5]) and using the conversation law of 
the total energy density p t + 3H(p + p) = 0, we can get 



H=--{p + p), 



(0) 



in which p = pde is the pressure of dark energy (the pres- 
sure of dark matter is zero), and pde = Pde^de-, with 
the equation of state parameter lode- We just consider 
lode is a negative constant, then it is another parameter. 
To discuss the dynamical behavior of interacting dark 
en ergy in EC, we introduce the dimensionless variables 



Pde 



y = 



Pdm 



Q 

3H 3 ' 



(7) 



Although the density of dark energy can be a negative 
on in the f(R) theory [22} > we restrict that the density 
of dark energy and dark matter both should be positive 
or one of then is zero. So both of x and y should be 
non-negative. 

The Friedmann equation (|S|) and Eq. © can be rewrit- 
ten as 

x + y=l, (8) 
H 3 / lodex \ 3 

W = -2{ 1 + 7T7j = -2 (1 + w) - (9) 

The effective total equation of state parameter u) is 

pDEWDE UDEX 



LO 



Pde + Pdm x + y 



ljdex. 



(10) 



In EC, the dynamical behavior of dimensionless vari- 
able x, y are 



x' = z — 3ujdex + SujdeX 2 , 
V = -z + 3u DE xy. 



(11) 
(12) 



in which prime denotes the derivative with respect to 
the e- folding number N — In a (ao = 1). Just as J?} 



in which z* is the value of z in the fixed points. If x* — 
0, the second subequation of Eq. fllUl) implies that z* = 
0. This is a point for pure dark matter dominated [3], 
we will not consider this condition. In this paper, we 
will discuss two interacting model: Qi = SHCiPde + 
3HC 2 Pdm, and Q 2 = i^H^C p n DM . C u C 2 and 
C are some constants. For simplicity, we will consider 
C\ = C 2 = C and n = 2 for those two models. Those 
models have been studied in EC by many authors 0, 
0-[l3|- But for convenience to compare the dynamical 
behavior of interacting dark energy in LQC and the ones 
in EC, we review the dynamical behaviors of dark energy 
with different interacting term in EC in the next two 
subsections. 



A. Model I: Q = 3HC\p D E + SHC 2 p D M 

This interacting term has been discussed in many pa- 
pers 0, 0-[l2| . In this subsection, for more simplicity, we 
just consider C\ = C 2 = C, the dimensionless interacting 
term is z = 2>C(x + y). Eq. (fT5|) can be rewritten as 

- u DE x*(l - x*) + C = 0, (14) 

It is easy to get the fixed points (x*,y*) from Eq. fbl]) : 



Point Ai 



Point B\ 



1 + ./1-4-2- l_./l_4-^- 



l-Jl-4-^- 1 + ,/l - 4-5=— 

1 "BE V U DE 



U5) 



(16) 



In order to determine the stability properties of the 
fixed points, it is necessary to expand the autonomous 
system X' = /(X) around the fixed points X», setting 
X = X* + U with U the perturbations of the variables 
considered as a column vector. Thus, for each fixed point, 
one can expand the equation for the perturbations up to 
the first order as U' = M • U, where the matrix M con- 
tains the coefficients of the perturbation equations. The 
stability property for each fixed point is determined by 
the eigenvalues of M [3] ■ According to this property, it 
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is easy to get the eigenvalues for the autonomous system: 



Point A\ 



Point _Bi 



Ai = -\/lo D e{lo D e - 4C) + -ujde, 



A2 = 3\/u D e(ojde - 4C), 



(17) 



Ai = --\Jlo D e{lode - 4C) + -w_d_b, 



A 2 = -2i\Jlode{ude - 4C). 



(18) 



The stability property for each fixed point is determined 
by the sign of Ai, A2. The fixed point is a stable if both 
of Ai, A2 are negative, and unstable if both of them are 
negative, if Ai, A2 have not the same sign, then the fixed 
point is a saddle point. To ensure that < < 1, 

< 4-^— < 1 should be satisfied. For Point Ai, it is 
always a saddle point for Ai < and A2 > 0. For Point 
Bi, it is a stable point always satisfied, for both of Ai, A2 
always are negative if < 4—2— < 1 is held. If total 

■' ° U1DE 

equation of state parameter ui = U)DE " PE v /l 4 "° e i s m 
the regions of— 1 < lo < — |, the universe will enter an 
accelerated stage, then Point Bi is an accelerated scaling 
attractor, this has been shown in the [||. If lo < — 1, it is 
easy to get H > from Eq.©, the universe will have a 
big rip. 



B. 



Model II: Q = 3 2_ " H 3 ~ 2rl C p n Dl 



This interacting term has been discussed by p|.[l4|. 
In this section, we just consider n = 2, the dimensionless 
variable is z = 3Cy 2 , then Eqs. (fT3l) can be rewritten as 

- lodex^I- x*) +C(l- x*f = 0. (19) 

The fixed points (x*,y*) for this model are 

Point A 2 : (1,0), (20) 

C LOde 



Point B2 



And the eigenvalues are 



LODE + C ' LODE + C 



(21) 



Point Ai 
Point B2 



Ai — 3u>de, A2 — 3u>de, (22) 
3lodeC 



Ai 



-ZlodEi A 2 



C + LODE 



(23) 



For Point A2 , the condition of x* , y» > and lode < 
constrain that C < and lode + C < should be held. 
It is always a stable point for lode < 0. The universe 
will entre a dark energy dominate stage, and lo = lode- 
For Point B2, C < should be satisfied, it is always a 
saddle point. So, for this interacting model, it has not 
any attractive behavior except the universe enter a dark 
energy dominated stage. 



III. DYNAMICAL BEHAVIOR IN LQC 

In this section, we discuss the dynamical behavior of 
interacting dark energy in LQC. As mentioned in the 



introduction, the Friedmann equation is modified in LQC 
and can be written as 



H 2 



1 



1 



P_ 

Pc 



(24) 



in which p c = 



Differentiating Eq. ([24"l) and us- 



«7 

ing the conversation law of the total energy density 
pt + 3H(p + p) =0, one can get 



H 



1 



(P + P) 1 



2p 

Pc 



(25) 



Considering the dimensionless variables ([7]), we can 
rewritten Eqs. (|24l25j) as 



(x + y) 
H 

IP ~ ' 



3# z 



(x + y) 



1 



xlode 



x 



y 



(2 



y)- 



(26) 
(27) 



The total effective equation of state parameter lo is 
Pde^de lodex 



Pde + Pdm 



y 



(28) 



Considering the dimensionless variables, the dynamical 
equations for our interesting system can be obtained: 



3x(l + lode) 
.3*11 + ^) (2-x 

XLODE 



V), (29) 
(2-z-y),(30) 



3y + 3y 1 



in which prime denotes the derivative with respect to the 
c-folding number N = In a (ao = 1). 

The fixed points (x*,y*) satisfy x' = 0, y' = 0. Simpli- 
fying Eqs. (|29|30|) . it is easy to get 



Case I 



Case II 



ir* + = 1, 
-?>lodex*{1 - 



+ z* = 0. 



(31a) 
(31b) 



a;» (1 + lode) + y* = 0, (32a) 
-3x*(l + lode) + z* =0. (32b) 



Comparing the above equations with Eqs. (ll3[) . we find 
that the expression for the first case is as same as in the 
EC. The Case II is a special case which just exists in LQC. 
This is because of the Friedmann equation has been mod- 
ified, then the dynamical equations in LQC are different 
from the ones in EC (see Eqs. (|29l30|) and Eqs. (|llll2p ). 
We will discuss these two cases for fixed point in the last 
section. 

As in EC, the values of fixed points depend on the 
special form of dimensionless interacting term z. We 
will consider two models: Q = ZHC\Pde + SHCiPdm, 
and Q = 3 2 ~ n H 3 ~ 2n Cp DM , with some constants 
C 1 ,C 2 ,C,n. 
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A. Model I: Q = ZHdpDE + 3HC 2 pdm 

In this model, the dimensionless variable z can be writ- 
ten as 



Q 

3H 3 



3dx + 3C 2 y, 



(33) 



in which C\, C 2 are some constants. 

Considering the above equation, Case I and Case II 
can be rewritten as 

uJdex\ - uJdex* + Cix* + C 2 (l - £*) = 0, (34) 
C 1 x* + C 2 [(-x*(1 + u; D e)] -x*(1 + w de ) =0(35) 

The first equation corresponds to Case I, and the second 
one corresponds to Case II. It is easy to find that x* = 
in the Case II. This means that y* = 0. The universe 
will enter other stuff dominated stage, this is beyond our 
interesting, so we will ignore this case in this subsection. 
As the last section, we just consider C\ = C 2 = C. 
Now, we consider Case I. Eg. (f3~4"|) becomes 



C 



ujde 

Then the fixed points (ic,,y*) are 



= 0. 



(36) 



1 



Point Ai, 



Point B3 : 



The eigenvalues are 



1 



1- a/1 



1-4- 



1-4- 



(37) 



(38) 



Point A3 



Point B-i 



Ai = -3 - -cjde - -v^deO^de 4 C)> 

(39) 



A 2 = 3^cj D e(w de - 4C), 



Ai = -3 - -CJDE + -V^DE^DE - 4C), 

(40) 



A 2 = -3^w DE (w DE - 4C). 



The density of dark energy and dark matter both are 
positive ones or one of them is zero, this condition given 
a constraint on the relationship between ujde and C. 
In this case, ujde < and < < \ both should 

be satisfied. For Point A3, it is a unstable point if 1 + 
woe + ^deC > and wde < — 2, and it is a saddle 
point if wde > —2 or 1 + ojde + ^deC* < and wde < 
—2. For Point B3, it is stable point if ujde > —2 and 
1 + lode + uJdeC > 0, and a saddle point if ujde > —2 
and 1 + ujde + ojdeC < or ujde < —2. 

Comparing the above fixed points and eigenvalues with 
the ones of classical scenario, we find the expressions for 
fixed points are same for it is not necessary to consider 
Case II, but the eigenvalues are different. The reason 



is that the eigenvalues are determined by the dynamical 
equations, the dynamical equations in LQC is different 
from the ones in EC. In EC, the singularity maybe hap- 
pens. But in LQC, for Friedmann equation has been 
modified, the singularity is replaced by the bounce al- 
ways. We will show this evidence in the next section. 



B. Model II: Q = S 2 -"CH 3 - 2n p^ M 



In the model, the dimensionless interacting term is z = 
3Cy n , then the Eqs. (|31l32j) can be rewritten as 

lodexI - ojdeX* + C(l - x*) n = 0, (41) 
(-l) n C(l +u D E) n x? - x*(l+u DE ) = 0. (42) 

To get the values of x„, we need the special value of n. 
As an example to use Case II to get the fixed points, we 
will consider n — 2. Then, the dimensionless variable z 
is z = 3Cy 2 . In this case, Eqs. (|41l42p can be rewritten 
as 

udexI — udex* + C(l — a;*) 2 = 0, (43a) 
Cxi (1 + u DE ) 2 - x, (1 + w DE ) = 0. (43b) 

Solving the above equations, it is easy to get the fixed 
points 



Point Ai : 
Point B4 : 

Point Ca : 



(1,0), 



c 



UJDE 



UJDE + C UJDE + C 
1 1 



c{i+uj DE y c 



(44) 
(45) 

(46) 



Noticed that — is also the solution of Eq. (l43bl) . 
but considering the first equation of Eq. (|32l) . should 
be zero when x* = 0, the universe will enter other stuff 
dominated, this is beyond our interesting. 
The eigenvalues for above points are 



Point A4 
Point B4 



Point C4 



Ai = -3(1 + ujde), A 2 = 
Ai = — Sujde, 

C + Cujde + UJDE 



3ujde, 



A 2 = -3 



UJDE + C 



(47) 
(48) 



Al — — -UJDE 

3 ^uj' de C 2 - 4C 2 - 4uj de C 2 - ACujde 
2 



C 



zUJde 



3 V^deC 2 ~ 4C 2 - AujdeC 2 - 4Cuj de 



C 



(49) 



For Point A 4 , the universe will enter a dark energy 
dominated stage. If —1 < ujde < e.g., the dark 
energy is quintessence-like type, this point is a stable 
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point, and it is an accelerated scaling attractor. If 
io DE < — 1) it is a saddle point. For Point B 4 , to en- 
sure 0<£*,2/*<l,C<0 should be satisfied. So it is a 
saddle point if C-\-Cude+uide < and a unstable point 
if C + Cio de + lo de > 0. For Point C 4 , C < 0,Wdb < -1 
and lo de C — AC — AujdeC — Alode < need to be satis- 
fied, it is a unstable point if C + Clo ee + lode < and 
a saddle point if C + Clode + lode > 0. 

So, this model in LQC will have an attractor behav- 
ior just when the universe enter a dark energy with 
— 1 < lode < — \ dominated stage. This is different 
from the one of EC, in which the dynamical system will 
have an attractor behavior when the universe is all types 
of dark energy dominated. As the first model, due to the 
quantum modification of Friedmann equation, the dy- 
namcial behavior for LQC is different from the one for 
EC. 



IV. NUMERICAL ANALYSIS 

In this section, we will analysis the dynamical behavior 
of interacting dark energy by using numerical tool. 

In this paper, we restricted that < < 1. This 

restriction gives a constraint on C,lode- As discussion 
before, although the expressions for the fixed points in 
LQC are as same as the ones in EC for Model I, but due to 
the different dynamical equations, the stability behaviors 
for each fixed points in LQC are different from the ones 
in EC. We show the stable regions in the parameter space 
{lode, C) for the Model I in Fig FT] It is easy to find that 
the stable region for LQC is smaller then the one for EC. 



c 




"V; 



FIG. 1: The stable regions in the (lode,C) parameter space 
for Model I. In EC, the stable regions are region II + III. In 
LQC, the stable region is region II. 

In Fig|2j we plot the evolution of the total equa- 
tion of state parameter Lo t (described by Eqs. (llOI28p ) for 
different models. As Eqs. (|10I28I) showed, the value of 
lo = PDE ^ DE is determined by the values ot Pde, Pdm, 
and also depends on the interacting term. But just as 
the diagram showed, the final value of lo will tend to 
constant. When lode = — 0.6, C = —0.1, we can get the 



total effective equation of state parameter — 1 < to < — ^ 
in the final data for Model I, not only in EC, but also in 
LQC. But for model II, lo will be less then —1, and the 
universe will be phantom-like dark energy dominated. In 
EC, the universe will have a singularity, but in LQC, the 
universe will enter an oscillatory region and will have not 
any singularity, just as Figl3] shows. 



\ (a) EC 


i 1 i 

- - - z=3C(x+y) 
z=3Cy 










i i i i 



2 4 6 8 10 

t 




2 4 6 8 10 

t 



FIG. 2: (Color online). The evolution to in LQC. The values 
of constants are C = —0.1, lode ~ —0.6, p c = 0.82. The upper 
panel is for EC, and the under one is for LQC. 

Figure 2] shows the stable point in the classical cosmol- 
ogy and in LQC for Model I. Although the expressions of 
the fixed points in the classical cosmology are as same as 
the ones in LQC, but due to the dynamical equations are 
different, so the stable region in (lode,G) space is very 
different, just as Fig.(l) shows, then the value of stable 
points for those two cosmology are different. 

In Fig03 we show the evolution of the densities in EC 
scenario and LQC scenario. In EC, the densities will 
become bigger and bigger. But in LQC scenario, the 
densities of dark energy and dark matter and the total 
matter will not divergence. The universe will entre a 
oscillating regime, and the singularity is avoided. 

V. DISCUSSION AND CONCLUSION 

In this paper, based on defining some dimensionless 
variables, we simplified the equations of fixed points. We 
found X* + y* = 1 for the fixed points should be sat- 
isfy both in EC and LQC. It is easy to understand this 
condition hold in EC, for x + y = 1 is nothing but the 
Friedmann equation, just as Eq.® shows. But in LQC, 
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(b) LQC 




) 1 


2 t 3 


4 




Pde 




- (b) LQC 


Pdm 






Pde + Pdm 

















30 
t 



FIG. 3: The evolution of Hubble parameter in the Model 
II. The values of constants are C = —0.05, lode — —0.6, p c = 
0.82. The upper panel is for EC, and the under one is for LQC. 




FIG. 4: (Color online). The phase diagrams of interact- 
ing dark energy. The upper panel is for EC, with lode = 
— 1.2, C = —0.1. And the under one is for LQC, with 
iode = -1, G = -0.1, p c = 0.82. 



FIG. 5: (Color online). The evolution of energy density p. 
The upper one is the evolutions of Pde,Pdm and the total 
energy density p = pr,E + Pdm in EC. The under one is for 
those variables in LQC. Both of those two diagrams are for 
LODE = -1.6, C= -0.1. 



if the quantum effect is dominated, it is impossible to be 
satisfied for p < p c . But notice that the quantum effect 
should be considered just when p is comparable to p c - 
So x* + y* = 1 means that our universe is described by 
general relativity very well when the dynamical system 
arrives its fixed point. The energy density is very far 
from the one should consider the quantum effect. It is 
easy to check that all fixed points of the interacting mod- 
els discussed in 11[ are determined by the condition of 
Case I. Case II is just held in LQC, but we can get this 
case from the conversation equation for energy density 



P : 



-3H{p DE (l + wde) + Pdm}- 



(50) 



When pde = Pdm = 0, this implies PdeQ + u DE ) + 
Pdm = and x' = y' = 0. Although the conversation for 
energy density should be satisfied both in EC and LQC, 
x^(l+u>DE)+y* = means that x^lojje = —1 in classical 
universe, but considering Eqs. (|llll2l) and setting x' = 
y' = 0, this is nothing but + y* = 1. So, Case II 
is just a special condition of Case I in the region where 
the universe is described by general relativity. When the 
universe is quantum effect dominated, x* + = 1 will 
not be satisfied, so, Case II is a special condition in LQC. 
Considering Eq. ([32|) . we can find that the fixed points are 
described by this case is happening when z* oc x™ with 
n 7^ 0, 1. Notice that, lude < —1 should be satisfied in 
this case for pde,Pdm > 0. Those fixed points which 
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are determined by Case II are staying in the quantum 
regions, and if they are stable points, the universe will 
not entre the Einstein cosmology. 

Due to the Fricdmann equation is modified in LQC, 
the dynamical equations in LQC arc different from the 
ones in EC. Under the help of the simplified equations 
of fixed points, we discussed the dynamical behavior for 
two different interacting models in LQC. We find that the 
stability of system dependent on the parameters ujde,C . 
Also, < < 1 gives a constraint on uirjE,C. The 

attractor solutions for autonomous system in LQC are 
very different from the ones in EC, this fact is obvious 
when we compare the stable regions in (ude,C) spaces 
of two conditions, although the attractor solution is ly- 
ing in the stage which is described by general relativity 



both for LQC and EC. The values of those parameters 
also should be constrained by the observation data, it is 
discuss by many authors [30l - l33l ] in EC, but it is still need 
more studying in LQC. Also, just as we expecting, the 
unstable region in LQC, the universe will entre an os- 
cillatory regime and avoid the singularity which always 
happens in EC. 
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